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We study and simulate jV — 2 supersymmetric Wess-Zumino models in one and two dimensions. 
For any choice of the lattice derivative, the theories can be made manifestly supersymmetric by 
adding appropriate improvement terms corresponding to discretizations of surface integrals. In 
particular, we check that fermionic and bosonic masses coincide and the unbroken Ward identi- 
ties are fulfilled to high accuracy. Equally good results for the effective masses can be obtained 
in a model with the SLAC derivative (even without improvement terms). In two dimensions we 
introduce a non-standard Wilson term in such a way that the discretization errors of the kinetic 
terms are only of order 0{a 2 ). Masses extracted from the corresponding manifestly supersym- 
metric model prove to approach their continuum values much quicker than those from a model 
containing the standard Wilson term. Again, a comparable enhancement can be achieved in a 
theory using the SLAC derivative. 
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1. Introduction 



Supersymmetry is nowadays an important ingredient in most theoretical developments of quantum 
field theory beyond the standard model. It allows for the unification of the three fundamental forces 
described by the standard model and is also incorporated in supergravity and string theory. In the 
low energy regime this symmetry is obviously not manifest and the question remains by which 
mechanism supersymmetry if realized in nature is broken. From non-renormalization theorems it 
is at least known that this has to be answered non-perturbatively. In this view the lattice might serve 
as an equally good approach as it has been before for gauge theories. However since supersymmetry 
is an extension of the Poincare symmetry of spacetime it is inherently broken on a spacetime lattice. 

Here we study and simulate jV = 2 Wess-Zumino models in one and two dimensions. Lattice 
theories with different lattice derivatives and discretization prescriptions which preserve parts of 
the supersymmetry are simulated. It is checked that fermionic and bosonic masses coincide and 
that unbroken Ward identities are fulfilled to high accuracy. By introducing a nonstandard Wilson 
term in the two-dimensional theory we can suppress common G{a) artifacts. To include dynamical 
fermions several algorithms are used and compared with each other. For a more thoroughfull 
presentation of our results we like to refer the reader to [[!]]. 



2. Quantum Mechanics 

In the continuum, the action of our first model is given by the action 

Scont = j dT^ 2 + ^W a + W+W''w) With W'(0) = ^M; (2.1) 

it is invariant under the following supersymmetric variations: 



8^^=£Y, 8Wy=-e(<j> + W'), 8^y = 0, 
= xj/e, 8Wy=(<j)-W')e, 8^\ff = 0. 



(2.2) 



In order to perform numerical simulations and compare with previously results [g] we have fix the 
potential to 

W(0) = f0 2 + ff- (2-3) 

A lattice version of this supersymmetric continuum theory raises a couple of questions. First 
we can ask whether the lattice model admits part of the continuum supersymmetry. Integrating out 
the fermions y and changing variables from the bosons <p to the so-called Nicolai variables 

£=<j)+W' (2.4) 

renders the bosonic continuum path integral purely Gaussian. Discretizing this sum of squares 

(2.5) 

— "naive 
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Figure 1: Left: The naive discretization of the continuum action fails to recover the correct supersymmettic 
continuum limit when Wilson fermions are used (green and red graphs). For every finite lattice spacing the 
extracted masses differ vastly from each other. The situation is different when SLAC fermions are used. Here 
the masses coincide within statistical bounds and tend towards the correct a — > limit. Right: The same 
graph as on the left, now using the improved actions as discussed in the main text. Here also the model with 
Wilson fermions exhibits degenerate masses at finite lattice spacing and tends towards the correct continuum 
limit. 



one easily verifies that one of the symmetries is preserved. Since the presence of the additional 
"surface" term improves the behavior of the action with respect to supersymmetry this action will 
be called Nicolai improved. 

Second we have investigated whether there is an optimal lattice prescription for the Dirac 
operator. In particular it is a well-known fact that (ultra-)local hermitean Dirac operators will 
introduce fermionic doublers thus spoiling the balance between bosonic and fermionic degrees of 
freedom. Two strategies might be pursued, to double the bosonic spectrum as well or to use the non- 
local SLAC derivative. The former requires then to amend the superpotential with a corresponding 
Wilson term while the latter is free of any such modifications. 

2.1 Degeneracy of mass spectra 

The most obvious physical consequence of supersymmetric theories is the degeneracy of masses 
between the bosonic and fermionic channels which is simply due to the fact that supersymmetry 
transforms corresponding states into each other. In Monte-Carlo simulations the masses of the 
lowest lying state can be read off from the exponential decay of the connected twopoint function 
which can be readily measured. For various lattice spacings a we have measured the masses for 
all models in both channels, see Fig. [l[ For all improved actions the presence of one unbroken 
supersymmetry suffices to find the degeneracy even at finite lattice spacing. However while naive 
Wilson fermions fail to recover the correct continuum limit as expected [^] and are still plagued by 
strong 0(a) artifacts for the improved action, SLAC fermions show considerably smaller deviations 
for finite a and are much less sensitive to improvement terms. 

2.2 Ward identities 

Another important check for the presence of supersymmetry in the lattice theory is given by the 
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Figure 2: Left: The Ward identities mentioned in the text as a function of jjc — yj. The green graph shows the 
"intact" identity while in red the deviation from the continuum result is clearly visible. The data were taken 
from a lattice with 20 lattice points using Wilson fermions. The analysis of SLAC fermions is somewhat 
more involved and all signals are worse by an order of magnitude than what is shown here. Right: The 
masses of the various two-dimensional models as a function of the lattice spacing. All spectra are degenerate 
in the scaling region, however lattice artifacts are much more pronounced for Wilson fermions than for SLAC 
or modified Wilson fermions. 

computation of several Ward identities. For any observable O and supersymmetry variation 8 one 
should find that 

8{O) = (8O} = (2.6) 
holds. With the particular choice for O = ^> x Wy an d 8 = 8^ we have checked explicitly the relation 

(VxYy) ~ (Uy) = (VxYy) ~ (Ufa + ^0) = 0, (2-7) 

the results are shown in green on the left of Fig. |[ On the other hand, since 8^ is not respected 
by the lattice action one might expect 

( Wx xj/ y ) - (M<!>y-W^) = (8®S) (2.8) 

to hold, cf. the left of Fig. @, too. 

By considering Ward identities we have seen that indeed one supersymmetry is preserved 
while the other is clearly broken at finite lattice spacing. Moreover the breaking of Ward Identities 
vanishes rapidly with decreasing lattice spacing and weak coupling g. 

3. Wess-Zumino model in two dimensions 

The action we start from now reads 

S C ont = J d 2 x(2d^d^ + ^\W'\ 2 + YMY), (3.1) 
where W' denotes the first derivative of the holomorphic superpotential and M is given by 

M = $ + W"P + + W"P-, P± = ^(l±y 3 ), V | . ; (3.2) 
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Again this model is invariant under the following set of supersymmetric variations: 

50 = \\r l Ei +£iy/ 1 , 50 = \j/ 2 e 2 + E 2 \j/ 2 
8\j/ 1 = _iw'£i + <90£ 2 , Sv/ 1 = -\W'£i -50£ 2 , (3-3) 
5y/ 2 = d§£i - \We 2 , 8y 2 = -<90£i - \We 2 - 

It is still possible to construct a local Nicolai map given by 

£=2(d$)+W',Z=2(d<l>)+W' (3.4) 
and the Nicolai improved bosonic action is thus 

5bos = ^(2(50),(^),+w;(^) x +w;(50). r +i|iv;i 2 ) 0.5) 

X 

while the fermionic part reads 

Sferm = ^MxyWy, M=M Q + W"(&)&xyP+ + W"(<fc)5^/>_. (3.6) 

x,y 

The chosen superpotential differs from the quantum mechanical one and now reads 

W(0) = ^0 2 + |0 3 . (3.7) 

This particular lattice actions leaves one of the four continuum supersymmetries intact. Unlike 
in the quantum mechanical case we consider only improved actions but choose different realizations 
of the Dirac operator. Standard Wilson fermions are certainly a natural choice since they are free 
of doublers and ultralocal and hence easy and fast to simulate. However, they suffer from large 
6(a) discretization errors and in our case necessitate a modification of the bosonic kinetic operator 
as well. SLAC fermions are again another choice. In order to obtain reasonable results, we have 
checked that the theory remains one-loop renormalizable via an explicit perturbative calculation. 
A third option emerges from a modification of the standard Wilson term reading 

(XT 

M = 7% + y /73A. (3.8) 

By this twist it can be shown for the free Dirac operator that all 6(d) artifacts vanish and the 
corrections become 6 '(a 2 ). Moreover on correlators of spatially averaged operators the corrections 
become even 6 (a 4 ) for the free theory. The right panel of Fig. || shows the masses of the lightest 
boson and fermion state respectively as a function of the lattice spacing. Both SLAC and twisted 
Wilson fermions are much less disturbed by lattice artifacts than standard Wilson fermions are, 
although the improvement of the action ensures the supersymmetric mass degeneracy in all three 
cases. 

4. Algorithms 

Since low-dimensional theories are less demanding than four-dimensional LQCD, several strategies 
to handle the fermion determinant on top of the standard hybrid Monte Carlo might be put to use. 
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Figure 3: Left: Distribution of the reweighting factor plotted as logarithm of the fermion determinant nor- 
malized to the free field determinant for different coupling strengths g using SLAC fermions on a 3 1 x 3 1 
lattice. The more pronounced the peak, the better statistical errors are under control and the more reliably 
estimates can be measured. It is obvious that the reweighting technique will fail for g > 1 . For each distribu- 
tion 20,000 configurations were evaluated. Right: Comparison of the bosonic two point function between 
the quenched and reweighted ensemble at g — 0.5 on a 32 x 32 lattice with Wilson fermions. The inclu- 
sion of fermionic fluctuations in the path integral are clearly vital for the correct computation of correlation 
functions and physical observables. 

In any case the models deviate at least in two points from the more familiar scenario of LQCD. 
First our theories involve a only single flavor in order to keep the fermionic an bosonic degrees of 
freedom balanced and second, y 5 -hermiticity is broken by the Yukawa coupling terms. In view of 
this and to start the investigation on safe grounds, various treatments of the fermion determinant are 
used in parallel and compared to each other. In the simplest case, the quantum-mechanical model 
with Wilson fermions, the explicit formula for the fermion determinant reads 

detM w [0]=n(l+m + 3£0 v 2 )-l, (4.1) 

X 

and can be applied directly to include fermionic contributions in a HMC integration scheme. Since 
the computational effort is rather small very high statistics are attainable. In the second related 
method one computes the determinant and the inverse of the fermion matrix by direct methods such 
as LU-factorization. Again additional noise originating from the use of pseudo-fermions is absent. 
While easily applicable in one dimension the method soon becomes infeasible in two dimensions as 
the lattices grow in size. A third possibility is given by reweighting the fermionic contribution from 
quenched ensembles. This method can generate configurations very quickly and is still exact in its 
treatment of fermionic fluctuations. Nonetheless it fails rapidly with increasing coupling constants 
since the fluctuations might then overstretch more than twenty orders of magnitude, see the left 
of Fig. ||, thereby reducing the effective number of configurations to order one. Finally pseudo- 
fermions are a well-known approach to estimate the fermion determinant stochastically. Recent 
algorithms such as PHMC and RHMC allow for the treatment of fractional powers of M^M. Thus 
these algorithms can also be used to simulate supersymmetric single-flavor field theories. However, 
the annoying problems with small eigenvalues of the fermion matrix will remain and may hamper 
the numerical treatment of these models. 
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5. Conclusions and Outlook 

We have tested several lattice constructions of supersymmetric = 2 Wess-Zumino models in 
one and two dimensions. The extended supersymmetry algebra admits the construction of a lattice 
action which preserves one supersymmetry. Using Wilson fermions this single remnant of the 
continuum symmetry suffices to observe important features of the theory such as a degenerate 
mass spectrum and the validity of associated Ward identities independent of the choosen lattice 
derivative. For the SLAC derivative in one dimension we have found that the results do not differ 
vastly between the naive and improved action respectively. With the help of the derived Ward 
identities it is possible to check explicitly that one supersymmetry is respected while the other is 
broken. 

The two-dimensional models are numerically more demanding since upon integrating out the 
fermion fields one ends up with an (in general not strictly positive) determinant. This situation 
worsens when the coupling is made stronger leaving this regime inaccessible for reweighting tech- 
niques. However the correct treatment of fermionic fluctuations is again crucial for the expected 
"supersymmetric" physics to show up as can be seen from the right of Fig. [| With the introduction 
of a modified Wilson term the typical 0{a) scaling is circumvented yielding results of about the 
same quality as the non-local SLAC fermions. 

In order to investigate the whole parameter space and/or models in more than two dimensions 
some technical obstacles related to the treatment of the fermion determinant must be readdressed. 
In particular we know from first experiments that preconditioning the linear systems before ap- 
plying iterative solver schemes would lead to a significant gain. Furthermore other acceleration 
techniques such as Fourier accelaration, multiple time-scales or higher order integrators are under 
investigation. With the help of the PHMC algorithm we hope to extend the stability of the algo- 
rithm into regions of parameter space which are inaccessible at the moment. Apart from this, it is 
already possible to study further supersymmetric models such as the j¥ = 1 Wess-Zumino model 
in two dimensions, nonlinear supersymmetric a-models, Wess-Zumino models in higher spacetime 
dimensions and Super- Yang-Mills theories with the help of our existing codes. 
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